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^ ■ Abstract 
Oh' 

D ' I review some aspects concerning the physics of neutrino mixing and oscillations. 

■ I discuss in some detail the physical neutrino oscillations parameter space in the 

\ case of two and three family mixing, and briefly describe the current knowledge 

of neutrino mixing parameters according to the present solar, atmospheric, and 
reactor neutrino data. I also briefly comment on the possibility of solving the LNSD 
, anomaly together with the solar and atmospheric ones. I conclude by emphasising 

\ that that even though in five to ten years time a lot will be learnt from the next 

round of neutrino experiments, a great deal about neutrino masses and neutrino 



. mixing will remain unknown. 
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1 Introduction and Motivation 



There are, currently, three neutrino puzzles, two which strongly indicate the 
presence of physics beyond the standard model (SM). The oldest one - the so- 
lar neutrino puzzle [1]- is the fact that the flux of electron-type solar neutrinos 
measured at the Earth is significantly smaller than predicted by solar physics 
models. This deficit is confirmed by different experiments, which make use of 
very different techniques for detecting neutrinos. It is fair to say that, after 
the publication of the first results from SNO [2], there is very strong evidence 
for a flux of v^^r coming from the Sun. 

The second - the atmospheric neutrino puzzle [3] - is the fact that the flux of 
atmospheric muon-type neutrinos and antineutrinos differs significantly from 
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theoretical predictions. This anomalous behaviour was first observed by the 
proton decay experiments 1MB and Kamiokande, and later confirmed and 
firmly established by the SuperKamiokande experiment. The anomaly mani- 
fests itself more strongly in the ratio of the predicted and flux-ratio to 
the experimentally measured one, and the nontrivial angular dependency of 
the flux. The latter is by far the most striking evidence for physics beyond 
the SM we have at the moment. 

Finally, the LSND experiment [4] , which studies neutrinos produced after pion 
and muon decays, has observes a 3-sigma excess of Pg-like events from /i+ 
decays. The LSND anomaly has not been confirmed or excluded by KARMEN 
[5] or other neutrino experiments. The situation will improve significantly with 
the advent of the MiniBoone [6] experiment, which is due to start taking data 
in 2002. 

The neutrino puzzles are best solved by assuming that the neutrinos have 
mass, and that neutrino mass eigenstates and weak eigenstates differ, hence 
the neutrinos "oscillate." The solar neutrino puzzle is best solved by assum- 
ing i/g <-> v^^r oscillations [7], while the atmospheric neutrino puzzle hints at 
quasi- maximal <-> i/r oscillations [8]. The LSND anomaly requires f>e 
oscillations. There are, however, other exotic solutions to individual neutrino 
puzzles in the market, such as new neutrino-matter interactions, violation of 
Lorentz invariance and quantum mechanics, and neutrino decays, to name a 
few. Usually, these solutions are tailor-made to address a particular neutrino 
puzzle. It is fair to say that neutrino oscillations is the only hypothesis that 
can properly address all neutrino anomalies. 

In this talk, I review some aspects of neutrino oscillations. In the next section, 
I discuss neutrino mixing and how it leads to neutrino flavour conversion, 
both in the absence and presence of a medium. I concentrate on two and 
three family mixing scenarios, and spend some time discussing the "physical 
parameter space." In Sec. 3, I briefly review how neutrino oscillations solve 
the neutrino puzzles, and what the current data can and cannot tell us about 
the oscillation parameters. In Sec. 4, I discuss some of the issues that must be 
faced if the LSND anomaly is also to be explained by neutrino oscillations, and 
in Sec. 5 1 conclude with an outlook regarding what we may expect to learn 
in the next few years and what will be left to do with future, more powerful, 
machines (neutrino factories?). 



2 Neutrino Mixing and Oscillations 

If neutrinos have mass, there is no reason for the mass eigenstates to coincide 
with the weak (also referred to as flavour) eigenstates. This not being the 
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case, the two bases are connected by a unitary matrix z/q, = UaiJ^i, where Ua 
are weak eigenstates (a = e, /i, r, . . .) and Ui are mass eigenstates, with masses 
rrii {i = 1,2,3, . . .). Uai will be referred to as the neutrino mixing matrix.^ 
Note that one can choose to label the mass eigenstates in ascending order of 
mass-squared, with no loss of generality. This will be assumed, unless otherwise 
noted. 

When neutrinos propagate in vacuum, it is very simple to calculate the proba- 
bility that a neutrino with energy E^, which is produced as a weak eigenstate 
Ua, is detected a distance L from the source as a weak eigenstate Ujs: 



where Am^^ = — mj are the mass-squared differences. For antineutrinos, 
one simply has to replace Uai ^ U*^^, meaning that eq. (1) holds as long as the 
sign of the last term is flipped. One can define the different oscillation lengths 



Figure 1 depicts L^{^ as a function of ^mf - for different values of E^. When 
the neutrinos propagate in matter (of either constant or varying density) the 
situation is significantly altered, as will be discussed in the next two subsec- 
tions. 



2.1 Two Neutrino Mixing 



If there are oscillation between only two neutrino states, the situation is rather 
simple. The neutrino mixing matrix is parametrised in terms of only one angle 
9, and one extra phase if the neutrinos are Majorana particles. This phase, 
however, is not observable in oscillation phenomena, and may be safely set to 
zero. Explicitly, if there is mixing only between z/g and z/^,, 

( ^ f cos 6* sin6'\/z/i\ 
\uj~\-sm9 cos9J\u2j' 

^ There is some controversy regarding how this matrix should be called. I will stir 
clear of this issue here. 
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Fig. 1. Lose as a function of 
Amfj, for different values of 
the neutrino energy, namely, 
E^=862 keV (^Be solar neu- 
trino), E^=W MeV ("typical" 
solar neutrino), Ey=l GeV 
( "typical" atmospheric neu- 
trino) and Ey=lQ GeV. 
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In the presence of neutrino-medium interactions, the neutrino propagation is 
modified [9]. These effects can be expressed in terms of a "matter potential," 
such that neutrino propagation is governed by the following Schodinger-like 
equation (in the weak basis). 



dL V^': 



2E„ 



sin^ 6 cos 6 sin 6 
cos 9 sin 9 cos^ 9 



+ 



V* V 



,(5) 



where V^fj are the different matter potential functions, which are in principle 
a function of the position. If is a linear combination of the other active 
neutrinos (z/^ and the matter potential matrix is such that Vex = and 
Vee — Vxx = V^GpNeJi^ if the neutrinos propagate in "normal" matter with 
electron number density A^e, and only standard model (SM) interactions are 
taken into account. For antineutrinos, is replaced by —Ng (which is the 
positron number density). 

Two important issues can be addressed just by inspecting the general form 
of eq. (5). First of all, while eq. (5) is, in principle, not solvable, we can still 
construct useful relations between the different Paf3- In particular, because the 

^ Note that the "Hamiltonian" in eq. (5) is defined up to a term proportional to 
the identity matrix, which does not yield any physically observable effect. 
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"Hamiltonian" is hermitian (this condition would break down, for example, in 
the presence of neutrino absorption by the medium or if the neutrinos decay), 
-Pa/3 = J2pPai3 = 1- This implies, in the case of two neutrino oscillations, 
that only one Pa/3 is independent. For example, one can express all Pa/s in 
terms of Pee' Pex — Pxe — ^ — Pee, Pxx — Pee- K is curious to note that, in the 
case of two family oscillation, Pex = Pxe (T-invariancc) is guaranteed, while 
nothing, in principle, prevents Pee 7^ Pee (matter induced CP- violation) . 

Second, it is interesting to define what is the "physical range" for the oscillation 
parameter 6, i.e. what are the values of 6 which yield, if one measures Pa/3{L), 
distinct results? Initially, because 6 is an angle, it need be defined only from 
to 2tt. On top of this, it is trivial to check that the transformation ^ tt — is 
equivalent to 6* — > —9, in the sense that the differential equations eq. (5) trans- 
form in exactly the same way. This implies that choosing values of ^ = [0, tt] 
is enough (this may also be seem by noting that eq. (5) depends only on 29, 
up to irrelevant terms proportional to the unit matrix). Furthermore, if the 
matter potential matrix is diagonal, the transformation 9 ^ tt — 9 modifies the 
differential equations in such a way that {ve^L), i'x{L))T^_g = (i/e(I/), —i/x{L))q. 
Because Pafj^L) = |z/^(L)p assuming the boundary condition is a pure Va 
state at L = 0, ^ and Tr — 9 yield the same PapiL), and the physical parameter 
space can be safely restricted to [0, 7r/2]. This is not that case if the mat- 
ter potential has off-diagonal terms and/or if the "weak" eigenstates at the 
production point differ from the "weak" eigenstates at the detection point. 
This is not the case within the SM, but in extensions to it one may indeed be 
able to distinguish 9 from n — 9. Finally, in the absence of a matter potential 
(which is the case of vacuum oscillations), there is yet another "symmetry:" 
(/^(L), 7^;(L))jr/2-0 = ('^*(-^)i ~^x(-^))ej such that the entire physical param- 
eter space can be spanned by allowing 9 to vary from to 7r/4. This can be 
explicitly seen in eq. (4). 

The parameter space in the absence of physics beyond the SM is easy to 
interpret. In the region 9 = [0,7r/4] (the light side) z/g is "predominantly light" 
(z/i), while when 9 = [7r/4,7r/2] (the dark side) Ue is "predominantly heavy" 
(z/2). While pure vacuum oscillations cannot tell the dark from the light side, 
matter effects brek the degeneracy and allow one to explore the entire [0, 7r/2] 
range. 

In many cases of interest, eq. (5) can be solved, at least approximately. If 
the neutrinos propagate in a constant electron number density. Pee is easily 
calculable, and is given by eq. (4), where the mass-squared difference and the 
mixing angle are modified to "matter" parameters, given by 

cos 29 m = I , (6) 

^JAl^ + A^- 2AA12 cos 29 + 43?(T4x) A12 sin 2^ + 4|KxP 
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Af^ = V^?2 + - 2AAi2 cos 29 + 43f?(Kx) Ai2 sin 29 + 4|KxP, (7) 

where Aj^ = ^rn^j/ {2E„) and A = Vge — Krx- discussed in the previous 
paragraph, the matter angle and frequency depend only on 29, and, if Vex = 0, 
are invariant under 9 ^ tt — 9. 

Another very useful solution exists in the case of neutrinos propagating in 
a monotonically falling electron number density, as is the case of neutrinos 
produced in the Sun's core [10]: 

Pee = Pi cos^ 9 + P2 sin^ 9 - cos 29m^Pc{1 - Pc) sin 29 cos(Ai2L + (^m), (8) 

where P^ = 1 - P^ = 1/2 + 1/2(1 - 2Pe) cos2^m, Pc is the "level crossing 
probability," (pM is a constant matter- induced phase, and the matter angle is 
to be computed at the production point. If the neutrino propagation inside 
the Sun is adiabatic, Pc = 0, while in the case of an exponential electron 
number density profile Ne{x) = A^"o cxp(— x/rg), the approximate form P^ = 
(exp(— 7 sin^ ^) — exp(— 7))/(l — cxp(— 7)), where 7 = 27rroAi2, is valid in a 
large portion of the parameter space [11]. Recently, a lot of progress has been 
made concerning the understanding of solar neutrino oscillations, including 
studies of whether the expressions above are valid for 9 > 7r/2 [12,13] (the 
dark side of the parameter space), and how they should be understood and 
interpreted [14]. 



2.2 Three Neutrino Mixing 

The three by three neutrino mixing matrix which connects Ve, i^n, and Ur to the 
mass eigenstates 1^1,2,3 is parametrised by three angles and one complex phase, 
plus two Majorana phases which are not detectable via neutrino oscillation 
experiments and can be safely set to zero. It is "traditional" to define the 
mixing angles ^12,13,23 in the following way: 

tan^ 9i2 = tan^ ^23 = sin^ = pesl^ (9) 

while 

Q{U*2Ue3U^2U*-i) = sin 6*12 cos 6*12 sin 6*23 cos 6*23 sin 6*13 cos^ 6*13 sin 5 (10) 
defines the CP-odd phase 5. 
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In the case of pure vacuum oscillations, the survival probabilities are given by 
eq. (1). It should be readily noted that, unlike the two-family case, if 5 7^ or 
71, Pa/3 7^ Paf3, siguaUing CP-violation in the neutrino sector.^ 

If the neutrinos propagate inside a medium, the situation is (much) more 
complicated. In principle, one has to solve the Schrodinger-like equation (in 
the weak basis) 

Up, where Vap = Vg*„. (11) 



On the other hand, because the "Hamiltonian" in eq. (11) is hermitian, the dif- 
ferent oscillation probabihties are related. As before, the constraints J2a Pap = 
1 and Y.pPap = 1 arise from the unitary evolution of the quantum state. In 
the case of three family oscillation, they imply that there are only four inde- 
pendent Pap, which one may choose to be, say, the three "diagonal" Paa plus 
Pe^- The other five are linear combinations of these four. Explicitly, 

P =1-P -P P =P +P -P 

^ er ^ ^ ee ^ e^i re ^ ~ e^ ^ m 

P^e 1 ~l~ Ptt Pee P^J^^^ -^e/i) (-^2) 

P =P +P -P P =l-P -P 

fiT ^ ee ' ^ e^ ^ tti ^ r/i ^ ^ ^ e^- 

Note that, unlike the two-family case, Pap may differ from Ppa-, even if the 
CP-odd phase 5 is "turned off" [15]. This means that the presence of the 
medium can induce not only matter CP-violation, but all matter T-violation. 
In order for this to happen, all that is required is that some Vap breaks T- 
invariance. This does not happen in the presence of a constant matter density, 
but may certainly happen, for example, for neutrinos produced inside the Sun 
[15]. Whether or not such effects are relevant for terrestrial neutrino beams 
has also been studied [16]. 

Also from eq. (11), it is possible to determine what is the "physical range" 
for the parameters 6ij and 5. One can directly check, as discussed in the two- 
family case, that if V^p = for a 7^ /5, V^^ = Vrr, and the weak eigenstates 
are the same at the production and detection locations (as are the conditions 
in "normal" matter, including no interactions beyond the SM), all physically 
distinguishable values for the oscillation probabilities are probed if the three 
angles are allowed to vary in the range [0,7r/2], while 6 covers the full range 
[— 7r,7r]. Note that, unlike the two-flavour case, there is no further reduction 
of the parameter space if the matter potential vanishes. In the presence of 
new physics effects that lead to a non-diagonal matter potential, a larger 

^ Of course, there is also T-violation, P^p 7^ Ppa, such that Pap = Pp^ (CPT). 
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parameter space is, in principle, required in order to describe all possible 
oscillation scenarios. 

One final convention issue should be mentioned. Thanks to the current ex- 
perimental data, it has become customary to redefine the order of the mass 
eigenstates if Am^2 > ^^23- Whenever this (an "inverted hierarchy") hap- 
pens, the mass-eigenstates are relabelled 1 — 3 — >• 2, such that Amfg 23 < 
and Am^2 < l^'^isl) always. This redefinition is done such that one can re- 
late 9i2 with the "solar" angle, 623 with the "atmospheric" angle, and 6*13 with 
the "reactor" angle, as will be discussed briefiy in the next section. Note that 
within this definition of the mass eigenstates {i.e., ml < m\ and is de- 
fined such that Am^2 < l^^ial)' ^^e sign of Am^a becomes "physical," and 
determines whether the neutrino masses are hierarchical, or whether 'm\,m\ 
are quasi-degenerate (Am^2 ^2) inverted hierarchy. It should be 

noted that oscillation experiments cannot distinguish whether all three mass 
eigenstates are quasi-degenerate or not. 

Since we are now dealing with three-by-three matrices, simple exact solu- 
tions to eq. (11) do not exist, even in the case of a constant matter potential 
(see, however, [17]). However, Nature seems to have been kind, and a few 
approximations seem to be well justified, including Am^2/I^"^23l lEI ^"^^ 
sin^ ^13 -C 1. Under these circumstances many of the well known two neutrino 
results can be applied. A good rule of thumb is that, in the presence of "nor- 
mal" matter (no new physics), ^13 is modified to a matter given by eq. (6) 
with 9 and A12 A13; 612 is also modified to a matter angle as defined 

in eq. (6) with 9 — > 9i2 and a modified matter potential: A — > A x cos^ ^13; ^23 
and 5 are not modified. For solar electron- type neutrinos (with Ey < 10 MeV), 
for example, 

= cos^ 9^3P!:{e ^9,2^A^A cos2 ^13) + sin^ ^13, (13) 

where P'^^ is given by eq. (8), and A13 effects are "averaged out" (see [9] and 
references therein). 



3 Neutrino Oscillations in Action 

Here, I briefly review what the experimental data has to say about the neutrino 
oscillation parameters if the neutrino puzzles are interpreted as evidence for 
neutrino oscillations. As mentioned in the introduction, I will ignore the yet 
unconflrmed LSND anomaly, and concentrate on addressing only the solar 

^ see talk by Eligio Lisi [8] concerning the violation of this condition. 
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and atmospheric neutrino data, while satisfying the constraints imposed by 
reactor neutrino experiments. 

First of all, reactor neutrino experiments [18] have, so far, failed to observe a 
depletion of the u^-^ux produced by nuclear reactors. The survival probability 



1 - sin^ 2^13 sin^ 



A23L 



(14) 



and, in hght of the atmospheric and solar data, Pgg ~ 1 implies sin^ 29is ^ 
0.1 and Am^2 ^ 10~^ eV^ (reactor neutrinos have energies of several MeV, 
while the most recent experiments probed L of order 1 km). Second, the solar 
neutrino puzzle requires values of Pge significantly different from 1, such that 
values of sin^ ^13 significantly different from 1 are required (see eq. (13)). This, 
combined with the reactor bound, implies sin^ ^13 < 0.1. 

Because the atmospheric data requires |Am23| ~ few x 10^'^ eV^, Am23- 
effects for solar neutrinos "average out", and the solar neutrino oscillations 
are sensitive to mostly Am^2 6^12 (hence these are referred to as the solar 
parameters). Currently, disjoined regions of the parameter space satisfy the 
data quite well [7], and they are referred to as: the large mixing angle region, 
which contains 0.2 < tan^ ^12 < 3, 10"^ eV^ < Amj^ < lO'^ eV^, the low- 
mass-squared-just-so region, which contains 0.1 < tan^ ^^12 ^ 10, 10^^° eV^ < 
Am^2 ~ 10~^ eV^, and the small mixing angle region, currently disfavoured 
with respect to the other two after the publication of the SNO [2,19,20] data, 
which contains 10"^ < tan^ ^12 < lO'^, 10"^ eV^ < Amj^ < 10"^ eV^. These 
regions are depicted in fig. 3 [20]. 

Finally, the atmospheric data requires P^^ close to 0.5 for "large enough" 
values of L/E^, while Pee is roughly 1. The latter is easily satisfied given the 
reactor constraints on Am^2 ^i^^ ^13 ) while the former determines the 
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values of jAmlg] and 623. Explicitly (ignoring matter effects), 

P^^ - 1 - (sin^ 2^23 cos" ^13 + sin^ 623 sin^ 2^13) sin^ (^|^) + 0(Ai2L).(15) 

Note that "maximal mixing" corresponds to sin^ ^23 — 0.5/(1 — sin^ ^13). De- 
tailed analyses of the atmospheric data in terms of neutrino oscillations con- 
strain 0.3 < tan^ ^23 < 3 and 1 X lO'^ eV^ < |Ami3| ^ 6 x lO'^ eV^ [8]. 



4 Including LSND — Four Neutrino Schemes 

The LSND anomaly may be interpreted as evidence for conversion. 
If the conversion mechanism is neutrino oscillations, the LSND result requires 
Pfie ~ 10~^ for neutrino energies of several tens of MeV and a baseline of 
roughly 30 metres [4]. In order to obtain a small enough oscillation length, 
values of Am^ > 1 eV^, much larger than the values required in order to 
address the solar and atmospheric puzzles, are required. If all three neutrino 
puzzles are interpreted in terms of neutrino oscillations, three different Am^ 
values are required, which implies the existence of (at least) four neutrino 
mass eigenstates. One the other hand, very precise LEP measurements at the 
Z°-boson pole indicate that there are only three neutrino species that couple 
to the Z^-boson. One is forced to conclude, therefore, that the LSND anomaly, 
combined with the solar and atmospheric puzzles, hints at the existence of a 
fourth neutrino, which is a SM singlet (hence referred to as "sterile," Ug). 

There are two rather distinct ways of organising the four neutrino masses- 
squared in order to try to solve all the neutrino anomalies. One is the "3+1" 
scheme, which has four hierarchical masses-squared, ml < m\ < m| < m|, 
such that Amf2 is responsible for the solar anomaly, Am|3 ~ Amf3 for the 
atmospheric anomaly, and Am34 ~ Ato|4 ~ Am^4 for the LSND anomaly 
(there are other variations, such as an "inverted hierarchy" for the masses 
of ^^1,2,3)- In this scheme, the 1/4 state turns out to be predominantly Ug, and 
the solar and atmospheric puzzles are interpreted in terms of (predominantly) 
active neutrino oscillations. The LSND anomaly, on the other hand, requires 
a small i/g.^t component in u^. A reasonable combined fit to all three puzzles 
exists, and is quite robust. The biggest obstacle to the 3+1 schemes is to 
satisfy bounds from terrestrial i^^ and i^e oscillation experiments, which are 
sensitive to the large value of Am^ required to solve the LSND anomaly. For 
more details, see the talk by Orlando Peres [21]. 

Another possibility is to "pair up" the neutrino masses, such that m| ~ ml ^ 
ml^2- In this case, Am^2 is responsible for the solar anomaly, Am|4 for the 
atmospheric anomaly, and Amfg ~ Am^4 ~ Am23 — Am24 for the LSND 
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anomaly (one variation is Amf2 responsible for the atmospheric anomaly and 
Aml^ for the solar one). These are referred to as "2+2" schemes. In this 
case, there are virtually no constraints from terrestrial neutrino oscillation 
searches but, instead, one has to worry about properly solving the solar and 
atmospheric neutrino puzzles. 

The problem is easy to understand. It has to do with the fact that, in a 
2+2 scheme, either solar Vf. or atmospheric Vp, oscillations have to be into a 
predominantly sterile state. Explicitely, let us assume C/e3 — t^e4 — 0, and 
Vpx = — (note that these conditions are weakly violated if one is to 
solve the LSND anomaly), such that 

vi — cos -dve + sin i9(cos C,Vr + sin (^Vs) , (16) 

V2 — — sin i?i/e + cos 'dicos + sin (us), (17) 

1/3 — cos0z/^ + sin 0(— sin + cosC^^s), (18) 

1/4 — — sin (pi/p + cos 4>{— sin + cos C^s), (19) 

In this case, (0) would play the part of the solar (atmospheric) angle, while ( 
controls the "amount" of sterile neutrino in both the solar and the atmospheric 
sectors. It is easy to see that if C = 0, z/^ (z/g) oscillates into a pure Ug (^r) 
state, while the situation is reversed if ^ = 7r/2. 

The current atmospheric data strongly disfavours (at more than 99% confi- 
dence level [3]) C = 0, while the solar data disfavours ( = 7r/2 [1] (this can also 
be seen in fig. 3). Detailed quantitative analysis [7,8,22] seem to indicate that, 
if all atmospheric and solar data are combined, acceptable fits can be found 
(the current best fit point is close to sin^ C = 0.2 [22]). It is curious to note 
that, according to [22], the point where the sterile component is evenly shared 
among the solar and atmospheric pairs (C = 7r/4) is strongly disfavoured by 
the data. As more neutrino data accumulates, it is not impossible to imagine 
that 2+2 schemes will be severly cosntrained by the data (perhaps as much, 
or more, than the 3+1 schemes). 



5 Conclusions and Outlook 

It is easy to summarise the results of sec. 3: | Am^gl and 623 have been rather 
well measured, while Am^2 and 612 have also been measured, but rather poorly. 
Nothing is known about the value of 613 (except that it is relatively small), 
the value of the CP-odd phase S, or the neutrino mass-hierarchy (the sign of 
Am|3 in the traditional parametrisation) . 

In the near future, the situation is bound to improve significantly, especially 
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in the solar sector. The on-going solar experiments (in particular SNO) will 
not only establish once and for all that neutrino conversion is the mechanism 
behind the solar neutrino problem, but will also improve our knowledge of the 
solar parameters. The near future experiments KamLAND [23] and Borexino 
[24] have the abihty to go even further, by not only breaking the current de- 
generacy in the {Ami2 x ^i2)-planc but also determining the solar parameters 
with unprecedented precision (KamLAND in the case of the large mixing an- 
gle solution [23,25,26], Borexino in the case of the low-mass-squared-just-so 
solution [12,27]). Furthermore, it should also be noted that the KamLAND 
reactor experiment offers the first realistic opportunity to observe a "real" os- 
cillatory pattern, while the Borexino solar neutrino experiment may establish, 
similar to what SNO will eventually accomplish, whether there are neutrinos 
other than i/g coming from the Sun [28]. 

In the "atmospheric sector," the ongoing K2K experiment [29], the future 
MINOS experiment [30], and the CERN CNGS project [31] will confirm at- 
mospheric neutrino oscillations, discover r-appearance, and extend slightly the 
sensitivity to 6*13, on top of precisely measuring the atmospheric parameters 
jAmlgl and ^23- 

Finally, after the MiniBoone data [6] is collected and analysed, the issue of 
whether the LSND anomaly is indeed a consequence of new neutrino physics 
will be settled. We will then be able to decide whether there is indeed a 
light sterile neutrino, such that all three neutrino anomalies arc explained by 
neutrino oscillations. It should be noted, as briefly discussed in the sec. 4, that 
the current neutrino data already places a significant amount of "pressure" 
on sterile neutrino oscillations, and the situation is bound to change as new 
experimental data from solar, reactor and long-baseline experiments becomes 
available (for an example, see [26]). 

In five to ten years times, the neutrino mixing matrix will still be far from 
well known: most likely, one of the mixing angles will simply not have been 
observed (unless the long-baseline experiments "get lucky"), and whether or 
not there is a nontrivial CP-odd phase will remain unknown. Furthermore, 
we will be unable to say whether the neutrinos masses are hierarchical {ml < 
ml < ml), or whether the hierarchy is inverted (mg < mf ~ m|, such that 
m| — mf <^ m|). 

These challenges have to be tackled by a new generation of neutrino exper- 
iments, which requires cleaner and much more intense neutrino beams, of 
different flavours (if possible, as is the case of a neutrino factory). 

I conclude by mentioning that there are other ways of probing the neutrino 
mixing matrix and mass spectrum. In particular, the future observation of 
neutrinos from a nearby supernova will probably shed light on some of the 
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issues raised above [32] , while continuing searches for neutrinoless double beta 
decay (which are the best probes of whether the neutrinos are Majorana or 
Dirac fermions) and a "kinematic" neutrino mass in the tritium beta-decay 
spectrum also add valuable information [33]. 
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